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Recently, Kob et al. pQ reported non-monotonic temperature dependence of a dynamic 
length scale, £ dyn , in a model supercooled fluid. Specifically, they found that £ dyn peaks 
around a temperature at which deviations from mode-coupling-like fits to the relaxation 
time become visible (T = 6.0 for their model) and that there is a decrease of £ dyn below 
the mode-coupling crossover temperature (T c = 5.2 for their model). According to Ref. 
PQ, the maximum of £ dyn signals a profound change in particles' motions. Its presence 
allows to identify the mode-coupling crossover without relying on any fitting procedure. It 
is consistent with a non-monotonic temperature dependence of finite size effects [2] . 

To determine £ dyn , Kob et al. froze all particles in a semi- infinite space in an equilibrium 
configuration and analyzed the dynamics of the particles in the other half-space as a function 
of the distance from the boundary. This "point-to-set" method to determine a dynamic 
length scale differs from earlier calculations [31 H] where the spatial correlations of particles' 
dynamics were analyzed. By examining these correlations, we found that the more commonly 
used |3j dynamic correlation length £ 4 increases monotonically with decreasing temperature 
in the range of temperatures investigated by Kob et al. 

We simulated the binary harmonic sphere system of Kob et al. using system sizes ranging 
from N = 10, 000 to 100, 000 particles with up to 10 9 time steps for the 100, 000 particle 
system at T = 5. To determine £4, we followed Refs. [U |5j. We defined a microscopic 
overlap function w n (t) = Q[a — \r n (t) — r n (0)|] where r n (t) is the position of particle n at a 
time t, Q(x) is the Heaviside step function, and a = 0.3 a (er is the diameter of the smaller 
particle). We calculated S±(q;t) = N^ 1 (^2 nm w n (t)w m (t) exp[q ■ (r n (0) — r m (0)]Y which is 
the structure factor calculated using the initial positions of the particles that have moved 
less than a distance a over a time t. Long range correlations of these particles are revealed 
by an increase of S4 (<?;£) at small q. To facilitate the fitting procedure used to determine 
£ 4 (t), we independently calculated x 4 (t) = lim g ^. £4(9; t) utilizing the method proposed by 
Berthier et al. [5] and described in detail in Ref. [5J. 

In Fig. [l](a) we compare ^4(t q ) (circles) and Kob et al.'s £ dyn (squares). The relaxation 
time r Q is defined through iV -1 (^2 n w n (r a )) = e" 1 . In contrast to £ dyn , the dynamic corre- 
lation length ^(tq,) is monotonically increasing with decreasing temperature. This is also 
evident from a direct examination of S4(q;T a ), Fig. [jjb): it is readily apparent that ^4(r a ) 
at T = 5 (circles) cannot be smaller than at T = 8 (triangles), as was reported for £ d y n . We 
shall emphasize that the finite wavevector data in Fig. 1(b) are ensemble- independent. As 
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discussed in Ref. [6], the integrals of the four-point correlation functions, i.e. the q=0 values 
in Fig. 1(b), are ensemble-dependent and this dependence is removed by adding correction 
terms derived in Ref. [5]. 

We do, however, find evidence for a change of the particles' dynamics in the temperature 
range examined by Kob et al, albeit at a slightly higher temperature, T « 8, than T « 6 
suggested in Ref. pp. By examining the relationship between the self- diffusion coefficient, 
r a , ^(tq,), and X4,( T a) we conclude that there is a crossover which manifests itself, inter alia, 
in a change from r a ~ exp[/c 2 ^4] to r a ~ exp^^ 2 ], Fig. ljc). 

We should emphasize the difference between £ dyn and £4. The former length pertains to a 
non-linear response of a fluid's dynamics to the pinning of one half of the particles whereas 
the latter characterizes the size of dynamic heterogeneities in an unperturbed bulk fluid. 
While a length evaluated from a linear response function can be easily related to behavior in 
an unperturbed fluid, such a relation is much less clear for the non-linear response. Ref. [1] 
and the present Correspondence highlight the importance of: (1) establishing a connection 
between non-linear response and dynamics in an unperturbed bulk fluid, (2) determining 
which length scale is the most relevant one for the understanding of glassy dynamics in the 
bulk, and (3) extending simulations beyond the mode-coupling crossover. 
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FIG. 1: The temperature dependence of dynamic correlation lengths, the four-point 
structure factor, and the relationship between the length and the relaxation time, (a) 

A comparison of the temperature dependence of £4(1"^) (circles) calculated in this work and ^ n 
(squares) from Fig. 3 of Ref. Q]. (b) The four-point structure factor 54(g;r Q ) (symbols) and 
the Ornstein-Zernicke fits, X4( r «)/[1 + {£,{ T a)q) 2 ], used to determine ^4,(r a ) (solid lines), (c) The 
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relationship between r a and ^4(r Q ) showing a crossover from r Q ~ e k i&( T a) i Q Tq ^ e fc 2£ 4 ( T a)_ 
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